Many-body interaction in semiconductors probed with 2D Fourier spectroscopy 
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A particular difficulty in studying many-body interactions in a solid is the absence of an ex- 
perimental technique that can directly probe their key characteristics. We show that 2D Fourier 
spectroscopy provides an efficient tool for the measurement of critical parameters describing the 
effect of many-body interactions on the optical response of semiconductors. We develop the basic 
microscopic theory of 2D Fourier spectroscopy of semiconductors in the framework of the three- 
band model (heavy holes, light holes, and electrons). The theory includes many-body correlations 
nonperturbatively and can be generalized straightforwardly in order to describe 2D Fourier spec- 
tra obtained in atomic physics. We establish a relation between the 2D Fourier spectrum and the 
many-body correlations. It is shown, in particular, that 2D Fourier spectroscopy provides a principal 
possibility to establish experimentally the origin of the fast decay of the memory term describing the 
Coulomb interaction between heavy- and light-hole excitons. The theory is applied to an analysis of 
the available experimental data. Experiments providing more detailed information are suggested. 
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Understanding many-body interactions in solids is one 
of the key problems of modern solid state physics (see, 
for instance, the recent review Ref. [II). The long-range 
Coulomb interaction between electrons leads to complex 
dynamics of the excitations in semiconductors and plays 
the principal role in the nonlinear optical response. Dif- 
ferent experimental techniques have been developed for 
studying the effects of many-body interactions. One of 
the most popular experiments is based on the four-wave 
mixing (FWM)2. In these experiments the sample is il- 
luminated by rays characterized by (non-parallel) wave 
vectors ki, k 2 , and k 3 . The outgoing signal is detected 
in the direction that corresponds to the nonlinear cou- 
pling of the excitation pulses, say, — ki + k 2 + k 3 . The 
advantage of such measurements is that the many-body 
contribution to the signal is not blurred by the strong 
linear (single-particle) component. The sensitivity of the 
FWM spectrum to the details of the interaction between 
the excitons and other many-body excitations makes it 
an efficient tool for probing the many-body properties. 
However, the standard FWM experiment does not allow 
one to make a distinction between different contributions 
of many-body interactions and correlations to the shape 
of the resonance. As a result it is difficult to interpret a 
FWM spectrum and to extract specific characteristics of 
the many-body interactions and correlations. 

Recently the more flexible technique of two- 
dimensional Fourier spectroscopy^^ has been applied 
to studying the semiconductor properties^. The gen- 
eral scheme is similar to that of standard time-resolved 
FWM experiments with three pulses propagating along 
ki, k 2 , and k 3 launched at t = t%, ti, and £3, respec- 
tively. The difference from a standard FWM experiment 
is that measurements are performed not at a fixed or just 
a few values of the delay time r = mm(i2>i3) — ii> but 
rather for a dense series of values lying in some inter- 
val. Subsequently, the Fourier transforms are done with 



respect to the delay time as well as with respect to the 
signal time. These two Fourier transforms constitute the 
two-dimensional Fourier spectrum. It is important to 
emphasize that in the experiments reported in Refs. 00 
the difficult problem of measuring both the real and the 
imaginary parts of the signal was resolved. These ex- 
periments, therefore, provide the information about the 
phase acquired during the delay time. As will be shown 
below, this allows one to make a distinction between the 
diffraction on the gratings created by the heavy- and the 
light-hole excitons. Thus, 2D Fourier spectroscopy ex- 
periments give substantial insight into the details of the 
many-body correlations and provide vital information, 
which is barely accessible using the standard approach. 
Despite many advantages provided by the 2D Fourier 
spectroscopy, the application of this technique suffers 
from the lack of the understanding of the spectra from the 
microscopic standpoint. Usually, 2D Fourier spectra are 
described in the framework of the phenomenology of the 
nonlinear susceptibility, which hides the relation between 
features of the spectra and the microscopic characteris- 
tics of the system. 

Here we develop the basic microscopic theory of 2D 
Fourier spectroscopy of the semiconductors in the frame- 
work of the three-band model (heavy holes, light holes, 
and electrons) . We use the perturbational approach with 
respect to the excitation field. The interaction between 
the excitons, on the contrary, is taken into account ex- 
actly. We show that the 2D Fourier spectrum gives the 
unique opportunity to measure the key quantities de- 
scribing the exciton-exciton interaction. We use our the- 
ory to analyze the experimental results of Ref. [6. Our 
calculations produce relations between the spectral fea- 
tures of the 2D Fourier spectrum and the parameters 
characterizing the many-body interaction that allow us 
to suggest experiments that would provide more detailed 
information. Because of the generality of the approach, 
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our theory is able to describe also the 2D Fourier spectra 
of other physical systems, such as molecular nanostruc- 
tures. 

The basic idea behind the derivation of the equations of 
motion of the exciton polarization is that the states with 
the definite number of particles are the eigenstates of the 
Hamiltonian of the unperturbed (i.e. without external 
field) semiconductor. Initially, before the first pulse hits 
the sample, the system is assumed to be in the vacuum 
state |0) with full valence and empty conduction bands. 
The excitation pulses couple the eigenstates of the sys- 
tem populating, thereby, states with different number of 
particles. The excitation pulses create directly pjp , the 
linear response polarization. This excitation is described 
in the rotating-wave approximation by 

(^ + t W « + r«)pW = _^(o (t)j (l) 

where I = 1,2,3 enumerates the exciton pulses according 
to their time order. Here and below the Greek letters 



k, A, fx, v are multi-indices denoting the exciton state 
according to K = {n K ,a K }, where n K is the type of the 
exciton equal to h or to 1 for heavy-hole and light-hole 
excitons, respectively, and a K is the helicity of the state. 
In Eq. (fTJ) we have introduced T K as the decay rate of 
the exciton state k, and uj r as the detuning, which is the 
difference between the frequency of the rotating frame 
and the exciton frequency uj k . The external source is 
specified by the dipole moment d and by the component 
of the i-th pulse, Eal(t), with the helicity a K . The lin- 
ear responses Pk\ in turn, serve as sources of the third- 
order polarization. The dynamics of the third-order po- 
larization is conveniently written in terms of the operator 
Dfj, u = [Bfj,, [B V ,H]], where H is the Hamiltonian of the 
unperturbed semiconductor and the operator B^ creates 
an exciton in the state fj,. The dynamics of the polariza- 
tion corresponding to the FWM signal in the direction 
— k/ + k m + k n , where I, to, and n enumerate the excita- 
tion pulses, is governed by 8 - 
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The parameters f3 and the memory functions F de- 
scribe the effect of the exciton-exciton correlations. Ob- 
viously, interactions involving higher number of parti- 
cles do not contribute to the third-order polarization. 
These parameters are defined by (3** = (D K \B^Bl^ and 
F^{t) = (e iHT D KX e- iHT Dl v ). The last term in Eq. © 
accounts for the Pauli blocking with the phase-space fill- 
ing parameters^ Cfy . 

As follows from Eqs. |T]) and ^ we need to take into 
account only the excitons with helicity of unit magni- 
tude, i.e. the ones coupled to the one-photon states of 
electromagnetic field. Therefore, in what follows, we will 
consider only excitons with a K = ±1. 

We would like to emphasize that Eq. |(3J) exactly 
accounts for the Coulomb interaction between exci- 
tons. The generality of the basic ideas and non- 
restrictive assumptions make this equation virtually 
model-independent. However, in order to present the 
relation between the 2D Fourier spectrum and the pa- 
rameters characterizing the many-body interaction in the 
most transparent way, we restrict ourselves to the short- 
memory approximation. 8 The effect of non- locality of the 
memory function on the 2D Fourier spectrum will be in- 
vestigated elsewhere. The short-memory approximation 
corresponds to relatively short biexciton life time^ or fast 
decay of the memory function^^ In the framework of 



this approximation the last term in the r.h.s. of Eq. © 
is substituted by an instantaneous term, which we absorb 
into the modified /3-parameter defining 

0$=0$+7$- (3) 

In what follows we will refer to as the modified mean- 
field parameter, outlining that it is local in time. 

The exact form of 7^ depends on the origin of the fast 
decay of the kernel. In the simplest case of short life-time 
one has& 

7 kA = F kX (0) (4) 

lfiU 2(r Al + r„) fiu{ ' { ' 

In the case of the fast decay of the memory function^iiS 
the modification of the /3-parameter takes the form 

^ = \{ D ^ H-i^Y v ) D ^)- (5) 

We study the polarization mixing induced by the 
many-body interaction using the standard three band 
(electrons, heavy, and light holes) semiconductor model, 
where the exciton destruction operator has the form 
B^ = J dxdy B^(x,y) with 

B»(x,y) = ^(x - y)vl Al) (y)c aA2) (x). (6) 



3 



Here <^(x — y) is the exciton envelope wave function, v 
and c are the annihilation operators acting on the states 
in the valence band and the conduction band, respec- 
tively. The indices cr K (i) specify the spin of the electron 
in the valence (i = 1) and conduction (i = 2) bands. In 
terms of B^(x, y) one can explicitly write^ 

D^u = I dx x . . ■dy 2 B tJi (xi,yi)B„(x2,y2)U(x 1 ,y 1 ; x 2 ,y 2 ), 

(7) 

where U is the energy of the electrostatic interaction be- 
tween two excitons. 

From the definitions of f3£* and 7** one can derive 
the spin selection rules. It can be shown that the con- 
tribution of the /3-term reduces to cx — i/3^\P K \ 2 P K — 
if3^\Ph\ 2 Pn, where R denotes the exciton state "conju- 
gate" to k, with the conjugation understood according 
to the rule {h, ±1} <-» {1, =Fl}. The 7-term turns out 
to be less restrictive having the form oc 7* |P.\| 2 -F«- 



For both fi and 7 the reduction of the indices is per- 
formed according to the same rule (3^ = /3^(2 — S K \). 
The expressions for (3 and 7 are combined together to the 
modified mean-field term according to Eq. ((3]). In what 
follows we treat f3^ as phenomenological parameters. It 
suffices for our purposes since we are interested in a re- 
lation between the 2D Fourier spectral features and the 
microscopic characteristics rather than in first-principle 
calculations of 2D Fourier spectra. 

We consider the 2D Fourier spectrum obtained in the 
so-called rephasing scheme^, when the excitation pulse 
corresponding to the conjugated field arrives first, i.e. 
I = 1 in Eqs. ([T]) and @. Resolving these equations with 
respect to the FWM polarization P R and performing the 
Fourier transform with respect to both the signal time 
t and the delay time r = min(i2,^3) — £1, we obtain for 
P K (w,fi) the expression 



P K =d|d| 2 ^^— - 



(w - w K + ir K )(w -uj k + iT K + 2ir A )(fi + w A + iT\) 



■ lUw.fi), 



(8) 
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where II K (w, fi) is the Pauli blocking contribution to the 
2D spectrum 
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iTx) 
(9) 

In Equations © and © the function /«(T) = 
e -T[i(w K -w A )+r K +r A ] describes the dependence of the 
spectrum on the time separation of the second and the 
third pulses T = £3 — ^2, the frequencies w and fi cor- 
respond to the signal time and the delay time, respec- 



tively, and ET } = / dtEr'(t). Deriving Eqs. © and 
([9|) we have used the assumption that the pulses are short 
compared to the characteristic dynamical time scales de- 
termined by w K and r K . The function g\(T max ) = 1 — 
e "- max (£2+WA+ir A )^ w it n r max being the maximal reached 
value of the delay time, accounts for the finite range of 
the delay time used in the experiments and explains the 
wavy character of the spectrum along the vertical axis. 

As follows from Eqs. |(5J| and © the spectrum has res- 
onances in the 2D (w, fi)-plane at points with the coor- 
dinates (w^, — w„). It is seen that a particular exciton 
state k produces only one resonance along the w-axis. 
The resonances along the fi-axis are produced by the 
nonlinear coupling of the exciton state k with different 
exciton states. Such a separation between the interaction 
with heavy- and light-hole excitons is possible solely due 
to the structure of the 2D Fourier spectrum. 

Simple comparison of Eqs. © and © shows that the 



2D spectra of the FWM signal created by the Pauli 
blocking and the Coulomb interaction between the ex- 
citons have qualitatively different form. The resonances 
on the spectrum produced by the Pauli blocking have 
the Lorentz form along the vertical and horizontal axes. 
The reason is that the dependence of the signal on the 
signal and delay time is essentially the free evolution of 
the polarization created by a short pulse. This evolu- 
tion has the form of oscillations, which produce a simple 
pole after the Fourier transform. At the same time the 
resonances created by the Coulomb interaction fall off 
asymptotically as oc 1/fi and as oc 1/w 2 , along the fi- 
and w-axes, respectively. This is the direct consequence 
of the fact that the FWM polarization is continuously 
excited by the polarizations of the linear response. The 
dependence on the signal time is found as a convolution 
of the respective Green function with the source. After 
the Fourier transform with respect to signal time it yields 
the product of Fourier images of the Green function and 
the source. This results in the asymptotic form oc 1/w 2 
because of the harmonic time dependence of these func- 
tions. 

More detailed consideration of the spectrum depends 
on the specific experimental situation. Keeping in mind 
the analysis of the experiments reported in Refs. we 
make several simplifying assumptions. First of all we 
note that the experimental spectra corresponding to the 
rephasing scheme have clear elongation along the vertical 
axis. According to the discussion above this means that 
the main contribution to the FWM spectrum in these 
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experiments comes from the Coulomb interaction. We 
employ this observation by neglecting the Pauli blocking 
term H k (oj,Q) in Eq. (jHJ. Additionally, we assume that 
the helicities of the excitation pulses and the detected 
signal are not resolved. Thus, the signal is obtained a 
sum of contributions ijHJ). Also, we assume that the ba- 
sic exciton characteristics do not depend on helicity, so 
that we have only two sets of parameters corresponding 
to light-hole and heavy- hole excitons. As a result, one 
ends up with four resonances in the (u>, i7)-plane plane, 
situated near the vertices of a square. Finally, we take 
T = 0. 



For a qualitative analysis we consider the situation 
when all resonances contribute to the 2D Fourier spec- 
trum independently, so that near the point with coor- 
dinates (ojfj,,— lo u ) we need to keep only the respective 
resonances in Eq. J5J). We denote this resonant contribu- 
tion by As follows from Eq. ([8]) and the assumption 
about non-resolved signal helicity £^ is obtained by sum- 
mation over different helicities of the excitons with the 
frequencies and w v . Employing the assumption of in- 
dependence of the material parameters on helicities, we 
obtain 



J 



4" 



(cj - ojfj, + iFfj,)(u - lu^ + iTfj, + 2iT v )(fl + uo u + iT v ) " 



(10) 



r 



where 

o" t i,ov 

(11) 

We note that as follows from this expression A 1 ^ = A™ v , 
i.e. these parameters depend only on the type of the exci- 
tons because of the summation over helicities. However, 
in order to improve readability of formulas, we keep the 
general notations where it is necessary. In Eq. (fTTj) I is 
the thickness of the sample, n is the refractive index of 
the material, and c is the speed of light. The magni- 
tude of £"{ui,£L) has the simplest form. It is a product 
of pole functions. The real and imaginary parts of the 
signal demonstrate a more complex structure. For exam- 
ple, the real part of the resonance, i?^, near the point 
— w„) changes its sign along the curve R"(uj,Q) = 0, 
which in a vicinity of — ui v ) can be approximated by 
a straight line 



i^(w,n) = arg(^) + 



2(r M + r„) i 
r,(r M + 2r,)^) + f:^ + ^' 

(12) 

We apply the developed description of the 2D Fourier 
spectrum for the analysis of the experimental data ob- 
tained in Refs. @. We would like to start from noting 
that, as follows from Eq. (p~2|) . the slope of the zero line 
R\(oj, f2) = is the constant —4/3 « —1.33. The experi- 
mental value of the slope is found to be equal to —1.3. We 
would like to stress the universality of the slope and to 
emphasize that such a good agreement with the experi- 
ment is provided by the fact that the (wh, — Wh)-resonance 
is the strongest one and the slope of the zero line is 
weakly affected by other resonances. The slope, however, 
strongly depends on the dynamical model describing the 
exciton polarization. In particular, if one takes into ac- 
count the Pauli blocking then the slope would take value 
depending on the relation between the Coulomb interac- 



TABLE I: Fitted values of the parameters of the system stud- 
ied in Ref. 0. The values of the exciton frequencies are pro- 
vided in the stationary frame. 
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FIG. 1: The real part (a) and the imaginary part (b) of the 2D 
Fourier spectrum obtained using the parameters from Table 
[I] These figures should be compared with Figs. 7e and 7f of 
Ref.B 



tion and the Pauli blocking contributions into the FWM 
spectrum. The slope takes the value —4/3 in the limit 
when the Pauli blocking can be neglected. 

From the experimental data of Ref. [6( we find the pa- 
rameters A^ using the least-square method. The exciton 
parameters lo^ and Y ^ are tuned in order to minimize the 
deviation between the theory and the experiment. The 
results of the fit are shown in Table HI The 2D Fourier 
spectrum corresponding to the fitted parameters is de- 
picted in Fig. [TJ 

The difficulty of extracting more detailed information 
from the results of Ref. |9 is that the measurements were 
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performed using linearly polarized pulses and the linearly 
polarized signal was detected. Much more detailed infor- 
mation can be obtained if the signals with circular po- 
larization, i.e. with fixed helicity, are used. As follows 
from Eq. (fTTj) a particular choice of the helicities of the 
excitation pulses and the detected signal allows a direct 
access to particular /3's. For example, if the excitation 
pulses have helicity a = +1 and a signal with a = +1 
is detected, then each resonance in (u), f2)-plane will be 
determined by the specific It is interesting to 

note that this provides the principal possibility to obtain 
experimentally the information about the origin of the 
fast decay of the memory term describing the Coulomb 
interaction between the heavy- and light-hole excitons. 
Indeed, this particular choice of the helicities excludes 
the contribution from the interaction between {h, +1} 
and {1, —1} excitons, so that only the term oc 7^'^ will 
mostly contribute. If the memory decay is caused by the 
short biexciton life-time then, as follows from Eq. ([4]), 
this should result in imaginary f3 i ^ . This will not be 
the case if the respective matrix elements of F£„ (t) decay 
faster^ than l/(T M + r„). 

We would like to note that even more flexible access 
to different matrix elements of 0" is provided by the ob- 
servational scheme with the time separating second and 
third pulses equal to the delay time, T = r. In this 
scheme the different sequences of excitation pulses, which 
enter symmetrically Eq. ((5J), turn out to produce differ- 
ent resonances along J7-axis. As a result, the 2D Fourier 
spectrum obtained using this scheme has four resonances 
along the vertical axis separated by io\ — lo^. 

We conclude by considering briefly the effect of inho- 
mogeneous broadening. The broadening is taken into ac- 
count by averaging the spectrum with respect to a joint 
distribution of the exciton frequencies. First, we note 
the different effect of the inhomogeneous broadening on 
the diagonal and off-diagonal resonances. While averag- 
ing is performed the resonances situated near (u^, -w^,) 
move along the diagonal resulting in elongating reso- 
nances in this direction. At the same time the width 



of the resonance in the direction perpendicular to the 
diagonal does not change and is determined by the ho- 
mogeneous linewidth. As follows from Eq. (TIT)]) , the 
half-width in the perpendicular direction <5 M of the mag- 
nitude of the signal near (w^, — w M ) is found from the 
equation (5% + 2Tl) 2 (Sl + 18r£) = 144r£ and is equal 
to Sfj, » O.SSr^. Applying this relation to the spec- 
tra obtained in Ref. [6j we obtain w 1.2meV and 
Ti w 1.4meV, which are in a agreement with the values 
found in Table HI 

The effect of the inhomogeneous broadening on the 
off-diagonal resonances is determined by the relation be- 
tween the distributions of the frequencies of the heavy- 
hole and light-hole excitons. For not too high values of 
the broadening the value of the light-hole - heavy-hole 
splitting can be considered to be fixed. As a result, av- 
eraging leads to elongation of the off-diagonal resonance 
along respective line, whose slope depends on the ratio 
between the values of inhomogeneous broadenings of the 
light-hole and heavy-hole excitons. 

In summary, we have developed the basic microscopic 
theory of 2D Fourier spectroscopy of semiconductors. We 
have shown that the resonant peculiarities in 2D Fourier 
spectrum are directly related to respective microscopic 
quantities describing the exciton-exciton interaction. Be- 
cause of the two-dimensional structure of the spectrum 
the contributions from the Coulomb interaction between 
different excitons (heavy- and light-hole) turn out to be 
naturally separated. We demonstrate, in that way, that 
2D Fourier spectroscopy provides a unique opportunity 
to extract the information regarding many-body corre- 
lations in semiconductors from direct measurements. In 
particular, we show that it is possible to obtain experi- 
mentally the information regarding the origin of the fast 
decay of the memory term describing the Coulomb inter- 
action between the heavy-hole and light-hole excitons. 
We have given a simple application of the theory analyz- 
ing the experimental data reported in Ref. [6j. 
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